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sitofsis 


The present work is aimed to verify the dynamic- 
shear angle approach for the prediction of stability chart. 

A flexible tool holder is designed and fabricated. 
The flexibility of the tool holder is much higher in the 
horizontal direction as compared to the vertical direction 
so that the machine tool workpiece system is assumed to be 
a single degree of freedom system. Response locus of the 
tool holder is determined by exciting the tool holder in 
the horizontal direction and measuring its response in the 
same direction. 

Experimental stability chart is obtained by 
machining hot finished seamless M.S. tube on L .B .25 Lathe 
using throw-away-tip -type carbide tools. 

A computer program is developed to solve the . 
equationsof threshold of stability. Theoretical stability 
chart is obtained using the computed results. 

Experimental stability chart shows a close 
agreement with theoretical results above cutting velocity 
70 m.p.m. when the builtup edge is insignificant. 
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IIITIIGDUCTXQh 

Metal cutting process is often associated with 
vibrations between tool and wo rk niece. These vibrations 
cause poor surface finish, reduce tool-life, affect machining 
accuracy and reduce production rate. 

Three types of vibrations occur in machine tools: 

(a) Forced vibrations are caused by the unbalancing 
affects in the machine tool work piece system. 

(b) Free vibrations are caused due to the shock 
transmitted through the foundation from the 
neighbouring machine s . 

(c) Self excited vibrations are also called chatter 
vibrations. Chatter is caused by interaction 
of cutting process with th n dynamic behaviour 
of the machine tool structure . 

Forced and fre ? vibrations are well understood 
and can be eliminated easily. Chatter vibrations pose a 
greater problem. These vibrations are due to the dynamic 
instability of the machine tool svstem. 

In order to predict the dynamic stability region 
of machine tool system, one has to analyse : 
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(a) Dynamics of machine Tool Structure, and 
(t) Djrnajaic s of Cutting process. 

1*1 Dynamics of .1 -aching Too l Stru cture 

Dynamic behaviour of a structure is represented by 
harmonic response locus of tbn structure. The response of 
the machine tool structure in the direction of applied 
harmonic exciting force is called the direct response of 
the structure. Direct response per unit force is termed 
direct receotan.ee. 

The response of the structure in the direction 
perpendicular to the applied harmonic exciting force is 
termed the cross response of the structure. Cross response 
per unit force is called the cross receptance. 

The receptance of a structure can be represented 
by a complex number, (a+ib). (?<>f. 1,2). Where a represents 
the inphase component of receptance and b represents the 
out of phase component oi the receptance. If (a+ib) is 
plotted in such a way that 1 a f is plotted as abcissa and 
'b' as ordinate, then the- harmonic receptance locus is 
obtained, 

1.2 Dynamics of Cutting Pro cess 

During chatter vibration, the tool vibrates 
relative to the work piece and cutting parameters such as 
depth of cut, rake angle and cutting speed etc*, undergo 



3 


cyclic vibration which lead to non-steady-stato metal 
cutting. 

T lusty and Polacek (?wi.3) assumed that variation 
in cutting forces depends only upon the chip thickness 
variation caused by relative vibration of tool and 
workpiece . 

Tobias and Fishnick (R~f .4) assumed that variation 
in cutting forces under dynamic cutting conditions depends 
upon chip thickness variation, cutting speed variation and 
tool penetration rate. 

Das (Ref. 5, 6) used the concept of universal 
machinability index for predicting the dynamic cutting 
forces, fie also considered that shear plane does not vary 
under dynamic cutting conditions. 

Knight (fief. 7) used the results of Das and Tobias 
(Ref .6) for the prediction of the stability of a simplified 
machine tool system. 

The urevious research workers have tried to 
predict the dynamic cutting forces from the characteristics 
of steady state cutting forces, however the theoretical 
analysis has not shown satisfactory agreement with expe- 
rimental results. 
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Sainth (Ref .8) gave a theoretical model for 
predicting dynamic forces from the- results of the steady 
state cutting tests. He considered the shear plane length 
as the most influential parameter in the steady state metal 
cutt ing. Considering dynamic shear angle, he gave a model 
(Fig.l) of dynamic cutting forces from which the expressions 
for incremental tangential cutting force dF^ and incremental 
thrust force dF t can he written as follows s 


dF = 1 Sin ( t + ’-f-‘ ) 
c c ' c 


( 1 . 1 ) 
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‘t ~ jjaplltude of the incremental thrust force 
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•‘ i c2 

h.± 2 

Sin 0 

(l - n s Cot 0) 

(1.11) 

/ , - 

— O 

f a — ri , /V / 


Y c2 

cu T 

1 - 2irs Cor, 0 /% 

(1.12) 

^2 = 

T q Vi ^2 
Sin 0 

s (l + 2c Cot 0) j 

(1.13) 


H 2 = r •Ws Cot Cf/* (i. 14) 


.'.vainth carried -out experiments to me asure dynamic 
cutting forces and showed good correlation between ex tie rim 
and theoretical results. 


Xsl — Determin ation of Threshold or S tability 

Under tho presence of chatter vibration 5 cutting 

forces (F c + dF c ) and (F^. + dF^.) a.ct on the system as show! 
in Fig. 2. 


Knight (Ref .7) gave the following expression for the 
displacement of tool in x -direct ion : 

X(t) = - (g + ill) dF t - (p + iq) dF c (l.l5) 

where , 

(g + ih) = direct receptance of structure 
(p + iq) = cross receptance of structure 
■It the threshold of stability, vibration magnitude 
will remain constant irrespective of time. 
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Hence vibration can be written in the form, 


XU) = e luJt 


= ds_ 


ClrlQ 


- tool vibration 
dSg ' = L q e iio( 'fc-T) 

- wave generated on the work piece 

, - . i W (t-T) + T 

dsp - u e 

* o 


Cl -iv 


(i.ifj 


value of dSg ! at a mean angular lead of "V" 
ahead of tool point 'O' 

Kainth (Ref. 8; used his analysis to solve equation 
( 1.16) and gave the following conditions for the threshold of 
stability : 

r 
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In equation CH9) 1 a'r’ imaginary parts can be put 
to zero separately 5 

g(l~ ACos © : )~k /'-Sin O’) | 0 - 11 ^( 1+20 Cot 0)J 

h 60 s ( , .. ■} 

j ■ “v ■ fviO Cot 0)-K 3 (l+2 c Cot 0)< 

1 c i 0 — j 

Sin0 


0 = 


jp(l-/i Cos 0 ! )-q /<- Sin O r j (i-r^s Cot 0) 
- q ^ S ' Cot 0 - cj 

[h(l-yUCos 6«)+ C .XL Sin j ^c-n-jSd+20 Cot 0)j| 

I 
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/ , i 

+ Jq (l- a*-'Zos 6‘) + p /a Sin © r j (l-n^s Cot 0) 

_+ t- U1-K-) Cot 0 - c j 


:i 
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CHAPTER 2 

EXFBR IhMTdL INVEST I GATION. OF THE STABILITY C&lRT 

m 

2.1 Introduction 

Present work is aimed to verify the dynamic shear 
angle approach for the prediction of stability of machine tools 
from the results of steady state cutting tests. For analytical 
simplicity a flexible tool holder is designed such that is 
essentially a single degree of freedom system. Direct harmonic 
response locus of the tool holder is obtained by exciting the 
tool holder by an electromagnetic shaker. Stability chart is 
determined by carrying out cutting tests on H.M.T. L.B.-25 Lathe. 

2.2 Flexible Tool Holder 

A flexible tool holder is designed and fabricated as 
shown in Fig. 3. Tool is secured in a block which is mounted 
on four leaf springs. Each leaf spring is screvred with the 
base of the tool holde£ at two points. This ensures rigid 
fixing of the leaf springs with the base of the tool holder. 
Other ends of the leaf springs are screx-rcd with the block 
carrying the tool. 

Head of the tool holder is joined to the block with 
the help of four bolts. By changing the mass of the head, 
supported on the leaf springs the natural frequency of the 
tool holder can be changed. 
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Leaf springs arc designed in a manner such that the 
section modulus of those springs in y-direction is sixteen 
times the section modulus in x -direction. This is done to 
ensure that the tool holder is essentially a single degree of 
freedom system. 


2.2.1 Response Locus of t he Tool holder 

Harmonic response locus of the tool holder is required 
for solving the equations of threshold of stability. This is 
determined by exciting the tool holder in x-dircction and mea- 
suring the output of the tool holder in the same direction. 

The block diagram of tho instrumentation used is shown 
in Fig. 4. Photograph of the set up is shown in Fig. 5. 

Tho electromagnetic shaker is placed in horizontal 
direction and is clamped to the base plato by a metal strip. 

To measure the force amplitude a dynamometer* is connected with 
the shaker through a right hand and left hand bolt. Lock nuts 
are provided to ensure the rigid fixing of dynamometer with the 
shaker, ii. mild steel bar threaded at both ends is used to 
connect the dynamometer with the tool holder. Tool holder is 
rigidly fixed with base plate through four bolts. 

Sinusoidal input to the shaker is provided by a 
function generator through a power amplifier. 


* Dynamometer designed by Loomba (Hef.9) was used. 
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L tri-axial accelerometer is fixed on the top of the 
tool holder to measure* vibrations of the tool holder in all 
the three directions. Output of the accelerometer is amplified 
in a charge amplifier and displayed on an oscilloscope. 

Dynamometer is calibrated statically by using standard 
weights. Calibration curve of the dynamometer is shown in 
Fig. 6. Cross response of the dynamometer in z direction 
was found to be negligible. 

Phase difference between the input and the output 
signals can be determined by a phase-meter provided both the 
signals are of required magnitude. The reference signal should 
be above 22 V, pk-pk. and the other signal termed as input 
signal should bo above 16 mV, pk-pk. In the present case 
dynamometer signal is of the order of a fraction of millivolts 
and accelerometer signal is of the order of a few hundred 
millivolts. Dynamometer signal is amplified by an amplifier 
having amplification factor of 400. The phase characteristic 
of the amplifier is determined seperately and is shown in 
Fig. 7. 

Phase difference between the dynamometer signal and 
the accelerometer signal is computed after measuring the phase 
difference of the two signals seperately with respect to a 
reference signal. Reference signal is taken from the function 

generator. 



msmmmm 
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2.2.2 Procedure 

111 the instruments are switched on, fer 10 to 15 
minutes, for stabilising the instruments b'forc conducting 
experiment , 

Force level of the shaker is kept constant at 7.5 kg.pk. 
which corresponds to 3 mV pk-pk output signal cf the dynamo- 
meter. The signal is displayed on the oscilloscope and its 


amplitude is adjusted by varying the input to the shaker. 

For a particular frequency, the force level is first 
adjusted. Output of the accelerometer in the direction of 
excitation is read from the oscilloscope. Inphase and 
quadrature components of the dynamometer signal and of the 
accelerometer signal with respect to the reference signal are 
measured scpera.tely by the phasemeter. The procedure is 
repeated by seeping the frequency in steps of 2.5 cps. 

Phase difference between the dynamometer signal and 


accelerometer signal is computed by substracting the measured 
phase differences of one signal from the other signal. Dis- 
placement. signal leads the acceleration signal by 180°. in 
account of it is made while computing the phase difference. 
Phase shift in the dynamometer amplifier is also taken into 


account. 

Ls the exciting sigral is harmonic, displacement is 
determined by dividing the amplitude of the accelerometer 
signal by 602 , where tu is the exciting frequency. 
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Tho output of tho accelerometer in the two directions 
other than tho direction cf excitation is below 5 percent of the 
output in the direction of excitation. Hones tho tool holder 
can be considered, essentially, a single degree cf freedom 
system. 

Inphaso component of tho accelerometer output is 
plotted as abcissa and tho out of phase component of the 
accelerometer signal is plotted as ordinate. The frequency of 
oscillation is indicated along tho curve. The curve thus 
obtained is called the harmonic response locus of the tool 
holder and is shown in Fig. 8. 

2.3 Experimental Stability Chart 

For determining stability chart experimentally, H.n.T. 
L.E.-25 Lathe being very sturdy is selected. Is the tool- 
holder is much more flexible in x-dircction as compared to 
y-direction, hence the machine tool work piece system can bo 
assumed to be a single degree of freedom system. 

Tho photographs of tho experimental set up arc shown 
in Fig. 9 and Fig. 10. 

For determining stability regions closer steps of 
speed, as compared to the steps available with L.B.-25, are 
desired. For obtaining closer stops of speeds, motor of the 
lathe is replaced by a continuous variable speed D.C .Motor. 








- SET UP FOB EXFEBIifflllTAL STiSILITY CSABT 
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A special mounting is designed said fabricated for this 
purpose . 


IxOt fiiiisnca seamless M.S. tube having .25% carbon 
was used for experiment. Carbide throw-away-tool-tips with 
5 rake angle wore used for machining tube specimens of 8' 1 


length at a feed of .05 mm per revolution. Higher feeds 
could not bo used due to the lower capacity of the variable 
snood motor as compared to the standard motor supplied with 
the lathe . 


2.3.1 Instrumentation 

Block diagram of the instrumentation used is shorn 
in Fig. 11. A vibration pick up is attached at the bad: of 
the tool holder. The pick up is suitable for a frequency 
range of 2-20000 cps . The pick up is connected to the 
vibration-meter for measuring the magnitude of the vibration. 
Output of the vibration motor is displayed on the storage 
oscilloscope for measuring its frequency. 

A techometor is used to measure the speed of the 
driving motor. Spindle speed of the lathe is computed by 
knowing the gear ratio and the speed of the original motor. 

2.3.2 Procedure 

Tube specimens of 8” length are turned with a no an 
radius of 54 mm. Work piece is mounted in the adopter which 
is used to avoid the centering of individual pieces. 






Thickness of tho tube is varied from ,5 mm to l mm in stops of 
.1 mm after oacn sot of tests for a range of cutting speeds 
at a particular thickness. 

ilmpl itude of vibration is measured by the help of 
tho vibration meter and tho frequency of vibration is obtained 
by displaying the signal on tho storage oscilloscope. Fresh 
cutting odge is used for each tost. After every reading 
the wavy surface left on t he workpiece, by the previous 
test, is removed by the help of back tool post. Three readings 
are taken for each speed for a particular width of cut. -Tests 
are conducted for spindle speeds of 100 to 500 rp.m. in steps 
of 20 rp.m. 

Figure 12 shows plots of freqaency versus revolutionary 
speed. Width of cut is also plotted against speed. The upper 
portion of the graph shows straight lines corresponding to 
stability lobes. Stability lobes are drawn touching tho points 
showing stability Such time all one points showing vibrations 
lie above the boundaries of the lobes. An envelope to those 
stability lobes is termed as threshold of stability curve. 












COPTER 3 





SfoLf Ifi tpo^uctlon 

In the present sthdy^ machine tool system is assumed 
a simplified single degree of freedom system* Hen.ce the 
cross receptance (p+iq) is taken zero. Also in orthogonal 
cutting value of overlap^factor ^ is one. Hence for the 


present case equationsof the threshold of stability reduce 
to the following equation.-, 

jg(l-Cos © l )-h Sin e'jjc-n^sCi+Rc Cot 0)j 

7 ((l+c Cot 0)-K-, ( l+2c Cot 0)? ra - 


I = - —S — - 
W Sin 0 


0 = 


I f If 

i jh (1- Cos ©•) + g Sin 9 T j jc-n-LS (l+2c Cot 0)j 


g^s 

V 


(1+ c Cot 0) rr. E 1 (1+ 2c Cot 0)j 


(3.2) 


whe re , 

©r = £sL (i 
n 


s Cot 0 
2 J\ v n ^ 


(3.3) 


Theoretical stability chart can be predicted by 
solving the equations (3.1) and (3.2). For this purpose, 
response locus of the tool-holder and the variation of the 
following cutting parameters with cutting speed is required : 
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(a) Cutting stress 'T f 

c 

(h) Shear angle , 0 r 

(c) Stress ratio 'D r and 

(d) Shear angle coefficients n x and rig. 

Variation of r T c ? , T D' 5 l n^* and 'ng* with cutting sreed 

Kainth (Ref .8) has given graphs for the variation 
of 1 T c 1 , *D ' , , n 1 ’ and 'Ug' with cutting speed. He also used 
hot finished M.S. tuhe with .25% carbon and carbide tool with 
5° rake angle. Cutting stress r T c ! , 'D 1 and T n 2 I are independent 
of chip thickness (Ref. 8; and their variation with cutting 
soeed is shown in Fig. 13 and 14. The variation of with 
cutting speed for chip thickness of ,05 mm is obtained by 
intrapolating the results of Kainth (Ref. 8). The intropolated 
curve is shown in Fig. 14. 

By using the multiple regression analysis, the 
following polynomials are fitted to the curves given in 
Figs. 13 and 14. 

T = 2.48 + 15,14 (V/10) - 1.1 (V/10) 2 + .04(V/l0) S 
c 

- .00043 (V/10)* 

D = .768 + 1.14 (V/10) - .78 ( V/100) 2 + .215 (V/100) 

- .02 (V/100) 4 
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n 1 = 1.348 (V/100)"* 2715 
n 0 = 3.4 (VX 3.281) "* 38 

o 

where T c = Cutting stress in kg/mm 2 
D = Stress Ratio 

n l = shear angle -chip thickness coefficient TraAiOWs/: ymrn. 
n 2 = Shear angle-rake angle coefficient 
V = cutting velocity in m.p.m. 

3.3 V ariation of Shear Angl e 0 wi th Cutting Speed V 

Experiments are carried out for determining the 
variation of shear angle with cutting speed. 

Mild steel specimen is machined to give a mean radius 
of 54 mm andtoidtfc of cut of 2.2 mm, Carbide tool with S b 
rake angle is used "or cutting, Feed of ,05 is kept 

constant . 

Method of length measurement is used to determine 
shear angle. A r Vee ( groove is cut on the tube, longitudinally. 
Chips collected during cutting carry no the s corresponding to 
tho groove. Chips equivalent of 4 to 5 revolutions are collec- 
ted and measured, thus giving the length of chip after machining. 
Tho notches on tho chips collected are used to determine the 
number of revolutions of the work-piece during cutting. 
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Therefore initial uncut chip length is, 

L =: 2 TT r N 

where r = Mean radius 'of work piece 

N = Number of revolutions for which chip length is 
measured 

If L 2 is the chip length measured, then the chip thickness 
ratio 'r c r is 



o s 2 I "| 

and shear angle 0 is, 

Xq Cos oi. 

0 = Arctan Sin , ;~) 

0 

where ^ is the rake angle of the tool. 

The variation of shear angle 0 with cutting speed 
is shown in Fig. 15. 

The following second order polynomial fs fitted to 
the above curve , 

0 = 6.482 + . S813(V/10) - .0067 (V/lO) 2 
whore 0 is in degrees and V is in m.p.m. 

3.4 Solution of t he S^-atlona-Qf Xte^hol^l^ 

For a particular frequency, roots of equation (3.2) 
are determined and for these roots values of widths of cut 
are calculated from equation (3,l)» 
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A computer program is developed to solve equations 
(3.1) and (3.2). The program is given in Appendix. 

So irch lor the roots of the equation 3.2 starts from 
60 r.p.iru for a particular frequency, 'Frequency range of 
160 c.p.s. to 230 c.p.s. is used in steps of 2.5 c.p.s. At 
th< st frequencies values of inphasc and out of phase components 
a and b respectively is determined from the hamonic response 
locus of the tool-holder. 

For a particular frequency function value of equation 
3.2 is calculated, at stops of 0,3 r.p.m. At each of the steps 
values of cutting parameters 0, D, n^ and n 2 are taken from 
the fitted curves. Whenever, at two consecutive steps, the 
function value changes its sign, a root is expected in the 

stop. Stocfanson's search technique is used to determine the 

-5 

root where the function value approaches 10 . All the roots 

are searched for a frequency and the whole procedure is repeated 
for all frequencies,, 

Values of cutting parameters T c , D, 0, n-^ and ng, 
at the roots determined from equation 3,2 is put in equation 
5.i and corresponding widths of cut are calculated. 

The theoretical stability chart is shorn in 

Fig. 16. 
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iAlSCUSSIOI i AND CONC LUSION 

A coranrrison of tho thoorotical and the experimental 
stability charts is shown in Figure 17. Experimental stability 
chart is a straight line in the speed range of 100 to 500 rpm. 
Vibrations of amplitude lessor than ,0125 mm arc considered 
stable. The tool holder is stable up to a width of cut .8 mm. 
For widths of cut of .9 mm and 1.0 mm., vibrations are of 
considerable amplitude. At those widths of cut certain 
speeds arc stable which lin in tho stability regions above the 
envelope of threshold of stability. Those stability regions' 
arc called regions of conditional stability. 

Theoretical stability chart shows a rising trend for 
speeds lower than 200 rpm. Between 220 rpm. and 500 rpm. the 
theoretical stability chart is quite close to the experimental 
stability chart, the maximum deviation Doing 12<5$>, 

At lower cutting speeds bolow 200 rpm, tho thoorotical 
values are higher than the experimental values. The reason 
is that Kainth's (Ref .8) theoretical analysis is not applicable 
for the built up edge conditions which is obtained below 
cutting velocity of 60 m.p.m. (180 r.p.m. in tno present cas^). 
However in practice carbide tools arc used at cutting velocity 
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hip'-nr tlnn 1,00 m.n.m. (300 r.p.nn), whore tho theoretical 

anal vs. Is shews f-'lrly rood agreement with tho experimental 
results. 

It Is conclud'd that tho dynamic shear angle 
a pn roach is ado quato to pr: diet the threshold of stability of 

machine tool system. 
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APPENDIX 


this program is to solve the equations of threshold q : stability 

$ I B JOB 
SIBFTC main 

REAL N1,N2,N,K1 

DIMENSION FF( 30 ) ,GG { 30 ) ,HH{ ?0) , ROOT ( 500) ,WW( 500 ) 
common /vishnu/theta ,ki ,c,v»fi ,s,ni , omega, h,f,gsmall 

COMMON /BLOCK /R00T»M 
500 FORMAT( 1X,10F10.3 ) 

510 FORMAT! IX , *FR EC'jE'J EMC Y* , 3 F 1 0 . n ) 

110 FORMAT < 8F10, 1 ) 

210 FORMAT ( 8F10. 6 ) 

R E AD 1 1 0 » { FF ( I ) ,1 = 1 ,30) 

READ?! 0 , (GG ( I ), 1 = 1 ,30) 

READ?l 0 , ( HH ( I ) , I = i ,30) 

DC inn 1=1,30 
F=FF ( I ) 

GSMAlL=GG( I ) 

H = HH( I ) 

X M I N = 2 • 0 
XMAX=10 .0 
5TEP=.005 
EPS=l.*E-06 

CALL.sTEEF ( XMI N , XM AX , STEP, EPS ) 

DO 10 J = 1 , M 
X=ROCT( J) 

CALL FUNCCX,Fl ) 

VI =V/ 1000. #60 • 

TC=10»# { •24'7S c 'iS + l • 5] 42A*V! /10.-. 10QR1 * ( Vl/10. )**?+. 003 58*<V]/10 
1 . )**?~.C0C04273*(Vl/10. ) **4 ) 

W = ~SIN ( FI } / T C / ( ( G s M A L L * ( l.-COS (THETA)) -H*SIN(THE tA))*( C-N1*S*Y 1 
1 . + 2.*C/TAN (FI ) ) )-H*O f 'EGA*S/V*( ( 1 ,+C/TAN ( F ID ) -Kl* ( l.+2«*C/TAN( FI ) ) ) ) 

WW ( J ) = W 

ROOT ( J ) =ROOT ( J ) *6 n , 

10 CONTINUE 

PRINT510,GSMALL,H,F 

PR I NT500 , ( ROOT ( J ) ,J = 1,M) 

PRINT500, ( WW ( J) ,J=1 , M ) 

100 CONTINUE 

STOP 
END 

SIBFTC FUNC 

SUBROUTINE FUNC(N,FF) 

REAL N 1 »N2 » Kl > N 

common /vishnu/theta, ki , c , v , f i , s , n i , omega, h,f , gsmall 

P I =4 , *A T AN ( 1 • 0 ) 

FD( X ) =.76R+1 •] 4’.63*X/100.-«7P.1R3* ( X / 1 00. ) **?_+• 21 55* ( X/100* ) 

1 **3-.0216R4*( X/10n. )*** 
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SI6FTC 

C 

C 

C 

c 

c 

c 

c 

c 


100 

101 

102 

103 


c 


10 

c 


FIA(X>=6.482 + .6812" 7 *(X/10. 3 — .006668* (X/10. 3**2 
FN1 <X)=1.348*(X/lft0# ) ** ( -.2715 3 
FN2(X)=3.4*(X*3.281 )**(-. 38) 

R = 54 . 0 
S = .05 

V=2 ,*P I *R*N 
V 1 = V/ 1 000 . *60 « 

M2=FN2 (VI) 

N1=FN1 (VI) 

FI=FlA(Vl)*PI/180. 

D=FD(V1 ) 

G=S/TAN (FI ) 

K1=G*N1+N2 

C=(D*COS(FI )-l.)/D/SIN(FI) 

OMEGA=2 .*PI*F 

T H E T A = OMEGA / N * ( 1 .-S/TAn (FD/2./PI/R) 

FF= ( { H* ( 1 • —COS C THETA) ) +GSMALL* S I N ( THETA ) } * ( C-N 1 *S* ( 1 . + 2 .*C/j AN ( F I ) 3 
) ) +GsmA|_L*OMEGA*S/V* ( ( 1 . + C/TAN (FI ) )-Kl*< 1 . + 2 .*C/TA,N ( f [ ) ) ) ) 

RETURN 

END 

STEEF 

SUBROUTINE STEEF ( XMIN,XMAX»STEP,EPS) 

PROGRAM IS For FINDING THE roots of TRAnsCIDEnTAl EqATIOM In the 


given range iteration function used requires no deravatio evalu 
notation used are 

XMlN=MlNlMUM value of InDErAnDEKT variable * 

xmax=maximum valuf of independent variable ***** 

STEP=STEP SIZE USED FOR SEARCHING THE VICINITY OF THE ROOT ***** 
ATION STILL CONVERGENCE IS faster THEN NEWTON'S METHOD ***** 

eps=accuracy desired For the function ***** 


DIMENSION ROOTS ( 5n0 ) .FUNCTN ( 500 ) 

COMMON/BLOCk/ROOTSjM 

F0RMAT(1CX»*N0. OF SOLUTION has EXCEEDED 20 *3 
FORMAT(10X»I4,2F12.5) 

FORMAT ( 10X,*SORRY NO SOLUTION IN THE GIVEN RANGE *3 

FORMAT ( 10X»*SOLUTION BY STEFFENSON'S IMPROVED REGULA r ALS I METHO 

D*//,10X»* N0.*5X ,*ROOTS' » 7X ,*FUNCT I ON*/ / ) 

SEARCH FOR THF VICINITY 6F THE ROOT **************************** 
X 1 =XM I N 

CALL FUNC(Xl.Fl) 

M = 0 

IF( ABS(Fl 3 .GT. EPS) GOTO 10 

IS SEARCH POINT VERY CLOSE TO THE ROOT ************************* 
X = X 1 
F = F1 

X2=X+5.*EPS 
CALL FUNC ( X 2 » F 2 ) 

G0T040 
CONTINUE 
X2=X1+STEP 
CALL FUNC ( X2 » F2 ) 

IS SEARCH POINT VERY CLOSE TO THE ROOT ************************* 
I F ( A B S ( F 2 3 .GT • EPS ) GOT020 




non 


F = F2 

X2=X+5.*EPS 
CALL FUNC(X2»F2) 

G0T040 

20 CONTINUE 

IF(Fi*F2.LT.0.0)GoTO30 

GOT050 

30 CONTINUE 

CALL PEGULA(X1,X2*F1,F2,EPS,X,F ) 

40 CONTINUE 

M=M+ 1 

IF(M.EQ.20)PRINT1o0 
ROOT S ( M ) =X 
FUNCtN(M)=F 

c if(m.eq .20 >dr inti pi , ( 1 .roots u ) , functn( 15,1 = 1 , 20 ) 

C IF (M. EG. 20 )M=0 

C SEARCH i- OR THE NEXT ROOT o T ART S ******************************* 

50 CONTINUE 

X1 = X? 

F1 = F2 

C TEST FOR THE SEARCH DOMAIN EXOUSTED OR NOT ********************* 

I F ( X 1 .LT.XMAX 5G0T01C 
IF (M.EQ.O) PRlNTlO? 

IF(M.GT.O) PRINT]03 

C IF^M.GT.O) PR I NT 101 , ( I , ROSTS ( I ) ,FUNCTN ( I ) , I =1 ,M) 

RETURN 

END 

$ I BFTC REGIJLA 

SUBROUTINE REGULA(X1 ,X2»F1,F2,EPS»X3,F3) 

100 FORMAT ( 10X »*NOT CONVERGING IN 30 ITERATION */) 

101 FORMAT ( 10X,*ITERATION NO.*, 13) 

102 FORMAT (10X,8F8.4) 

STEFFENSON'S IMPROVED REGULA FALSI ITERATION FUNCTION ***** 

THIS ITERATION FUNCTION IS OF THE ORDER OF 3*23 REQUIRS NO ***** 
DERIVATIVE EVALUATION. FUNCTION IS WITH MEMORY **************** 
M = 1 

10 CONTINUE 

X3=(X1*F2-X2*F1 ) / { F 2 — F 1 ) 

IF(F]*F2.LT.0.0)X3=(Xl*X3-X2*-*2)/(Xl-2.*X2 + X3) 

CALL FUNC(X3,F3) 

I F ( ABS ( F 3 ) .LT.EPS5GOTO60 
IF(M.EG.30 )G0T050 
20 CONTINUE 

C PRINT102,X1,F1,X2,F2,X3,F3 

IF(Fi*F3.LT.0.0)GOTO30 
X 1=X3 , 

F1=F3 

G0TO40 

30 CONTINUE 

F? = F3 


50 




X2 = X 3 

40 

CONTINUE 

M = M+1 

GOTO 10 

50 

CONTINUE 

C 

PRINT100 

60 

CONTJNUF 

X2=X3+5 , 0*EPS 
CALL' FUNC(X2,F2 

C 

PRINtIOI ,m 
RETURN 

SFNTRY 

13 

END 



